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SUMMARY 


Two  closely  related  numerical  methods  which  employ  operations  tables 
have  been  developed  and  used  in  the  calculation  of  the  buckling  load  of 
a  monocoque  cylinder  subjected  to  pure  bending.  They  are  based  on  the 
assumption  of  a  simplified  structure  which  includes  only  the  most  highly 
compressed  portion  of  the  cylinder.  The  first  method  makes  use  of  a 
Ih  -row  determinant j  whereas  the  second  method  requires  the  solution  of 
a  single  10-row  determinant.  The  buckling  loads  of  three  cylinders  with 
widely  different  characteristics  were  calculated  by  these  methods. 
Reasonable  agreement  with  experiment  was  obtained. 

A  procedure  similar  to  the  first  method  was  developed  for  the  cal¬ 
culation  of  the  buckling  load  of  a  cylinder  with  a  cutout.  A  limited 
experimental  check  was  obtained. 


INTRODUCTION 


The  calculation  of  the  buckling  loads  of  reinforced  monocoque 
cylinders  is  a  problem  of  some  importance  in  airplane  stress  analysis. 
Existing  theoretical  methods  for  detennining  such  buckling  loads, 

Including  energy  methods,  are,  in  general,  lengthy  and  difficult  to  apply. 

A  numerical  procedure  is  therefore  developed  in  this  report  in  order  to 
simplify  the  calculations . 

Southwell's  relaxation  procedure  (reference  1)  and,  in  general,  methods 
which  make  use  of  an  operations  table  (see  appendix  A)  have  been  success¬ 
ful  in  the  solution  of  a  variety,  of  stress-distribution  problems.  It  was 
therefore  natural  that  an  attempt  be  made  to  adapt  these  methods  to 
buckling-load  calculations.  In  reference  2  three  closely  related  methods 
for  determining  the  buckling  load  from  an  operations  table  were  established 
and  described.  A  limited  experimental  check  was  also  obtained.  In  refer¬ 
ence  2  the  three  methods  were  called  the  Determinant,  Energy,  and 
Convergence  Methods.  In  this  paper,  the  first  two  of  these  methods,  along 
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with  slight  modifications,  are  used  to  calculate  the  buckling  load^in 
pure  bending  of  four  monocoque  cylinders  with  widely  different  charac¬ 
teristics,  one  of  which  had  a  symmetric  cutout  on  the  compression  side. 

It  was  found  that  the  buckling  loads  could  be  conveniently  calculated 
when  the  actual  cylinder  was  replaced  by  a  simplified  structure  pre¬ 
serving  the  main  characteristics  of  the  original  cylinder. 

A  twofold  purpose  is  thus  fulfilled  by  this  investigation.  In  the 
first  place,  a  method  which  is  fairly  short  and  reasonably  accurate  is 
developed  for  the  calculation  of  the  buckling  load  of  a  monocoque  cylinder. 
Secondly,  a  further  experimental  check  of  the  methods  of  reference  2  is 
afforded  by  a  comparison  of  the  theoretical  and  experimental  buckling 
loads  for  the  cylinders  considered. 

The  authors  are  indebted  to  Dr.  N.  J.  Hoff  for  his  advice  and  help¬ 
ful  criticism,  and  to  Messrs.  J.  Mele,  B.  Erickson,  and  E.  B.  Beck  for 
their  part  in  the  experimental  phase  of  the  investigation. ^  The  work  was 
sponsored  by  and  conducted  with  financial  aid  from  the  National  Advisory 
Committee  for  Aeronautics. 


CALCULATION  OF  BUCKLING  LOAD  OF  CYLINDERS  WITHOUT  CUTOUT 
Methods  of  Calculation 


The  buckling  loads  were  calculated  for  three  cylinders,  the  charac  ^ 

teristics  of  which  are  given  in  table  I  and  figure  1.  The  methods  of 
calculation  which  appeared  most  convenient  are  described  below.  These 
methods  yield  the  load  P  in  the  most  highly  compressed  stringer  at  the 
instant  of  buckling  of  the  cylinder  as  a  whole.  From  this  load  the  total 
applied  bending  moment  can  be  calculated,  provided  the  stress  distribution 
is  known.  The  validity  and  the  accuracy  of  the  methods  are  discussed  in 
the  next  section.  Basic  theoretical  considerations  underlying  the  cal¬ 
culations  may  be  found  in  reference  2  and  in  appendix  A.  A  numerical 
example  is  given  in  appendix  B. 

Simplified-cylinder  solution.-  Let  the  cylinder  under  consideration 
be  replaced  by  the  simplified  structure  of  figure  2.  The  operations 
table  corresponding  to  this  structure  is  that  presented  in  table  II.  All 
symbols  which  appear  in  this  table  are  defined  in  appendix  C.  It  should 
be  noticed  that  the  operations  table  is  symmetrical  about  its  main  diagonal. 

The  buckling  load  P  has  the  value  which  will  make  the  determinant  repre¬ 
sented  by  table  II  equal  to  zero.  It  may  be  most  conveniently  obtained  by 
evaluating  numerically  the  determinant  for  several  values  of  P,  plotting  ^ 

the  determinant  values  against  P,  and  reading  off  the  load  at  which  the 
determinant  is  zero.  If  the  load  P  is  lower  than  the  first  buckling 
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load,  the  determinant  will  be  positive  (because  it  contains  an  even 
nmiber  of  rows;  see  reference  2).  If  the  determinants  are  evaluated  by 
the  method  of  reference  3,  that  portion  of  the  so-called  "auxiliary  • 
matrix"  which  corresponds  to  the  first  nine  rows  of  table  II  need  be 
considered  only  once,  since  it  is  independent  of  the  load  P. 

Solution  with  assumed  displacements.-  The  above  method  can  be  sim¬ 
plified  by  assuming  the  following  expressions  for  the  radial  displace¬ 
ments  r  and  the  rotations  of  the  most  highly  compressed  stringer 

(stringer  1,  fig.  2): 

r  =  sin^(nx/6L) 

m^  =  (dr/dx)  =  (5n/6L)  sin^(nx/6L)  cos  (nx/6L) 


in  which  the  maximum  radial  displacement  is  taken  as  unity,  and  L  is 
the  ring  spacing.  At  rings  B,  C,  and  D,  respectively  (see  fig.  2) 

+  determinant  is  then  reduced  to  that  given  in 

presentation  of  this  table  advantage  was  taken  of 
s^etry.  The  buckling  load  may  be  obtained  from  this  determinant  by 
e  first  method  given  above.  Since,  however,  only  the  element  in  the 

w^r+Iif^h  P,  the  determinant 

willptake  the  form  CK+f(kL)],  where  K  is  a  constant  not  dependent 


f(kL)  is  given  in  table  III,  and  (£1)3^^^  is  the  radial  bending 

stringer  and  its  effective  width  of  sheet.  The  value 
0  kL  at  buckling  makes  the  determinant  vanish  and  may  be  obtained 
irom  the  equation 

f(kL)  =  -K  (3) 

solw^thi°s  figure  3  and  may  be  used  to 

solve  this  equation  in  a  convenient  manner.  Consequently  the  buckling 

load  of  a  cylinder  can  be  obtained  from  the  solution  of  a  single  10-row 


the  and  lower  limit  may  be  found  for 

the  value  of  kL  at  buckling,  such  that 


l.ii6  <  kL  < 


h 
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The  value  1;.U9  corresponds  to  the  lowest  load  at  which  a  main-diagonal 
element  (in  the  tenth  or  eleventh  rows)  of  table  n  becomes  zero.  The 
value  1.U6  is  the  value  at  which  f(kL)  =0  and  is  approximate  since  it 
depends  on  the  assumption  of  equations  (1). 

It  mav  be  noticed  that  both  methods  require  the  evaluation  of  at 
least  one  determinant.  It  is  suggested  that  this  evaluation  be  carried 
out  by  the  method  of  reference  3-  The  following  remarks  concerning  the 
application  of  this  method  in  the  present  problem  may  be  useful: 

(1)  The  operations  table  is  symmetric  about  its  main  diagonal 

(2)  The  value  of  the  determinant  is  equal  to  the  product  of  the 

main-diagonal  elements  of  the  auxiliary  matrix  (defined  in 
reference  3) 

(3)  The  determinant  will  be  equal  to  zero  when  the  last  ^in- 

diagonal  element  of  the  auxiliary  matrix  vanishes  (see 
appendix  A) 


Discussion  of  Methods 

The  methods  outlined  in  the  preceding  section  are  based  on  the  ^ 
simplified  structure  of  figure  2.  The  following  considerations  underlie 
the  choice  of  this  structure  and  of  the  methods  of  calculation: 

(1)  The  most  highly  compressed  stringer  was  considered  of  paramount 
importance  at  buckling,  so  that  it  was  thought  permissible  to  neglect 

all  other  stringers  in  these  approximate  calculations.  This  is  equivalent 
to  considering  the  most  highly  compressed  stringer  as  a  column  elastically 
supported  by  the  rings  and  sheet.  The  elasticity _ of  the  supports  is 
represented  by  the  ring  and  sheet  influence  coefficients  in  the  operations 

table  (appendix  B). 

(2)  Points  on  the  tension  side  of  the  cylinder  will  undergo  o^y 
negligible  displacements  and  hence  may  be  considered  fixed.  The  rings 
are  therefore  assumed  to  continue  up  to  a  point,  near  the  tension  side, 

90°  away  from  the  most  highly  compressed  stringer,  and  to  be  rigidly 
fixed  there  (fig.  2). 

(3)  It  would  seem  natural  to  continue  the  sheet  up  to  the  same  point 
as  the  rings.  Because  all  stringers  except  the  most  highly  compressed 
one  have  been  neglected,  this  would  imply  a  single  panel  of  sheet  in 
each  bay,  extending  over  90O.  The  operations  table,  however  is  set  up 
considering  each  panel  with  its  edge  reinforcements  as  a  unit  in  which 
only  the  corner  points  have  independent  freedom  of  motion  (see,  e.g., 
reference  1|).  Therefore  the  action  of  the  90°  sheet  panel  would  be 
determined  by  the  displacements  of  its  corners,  with  no  possibility  of 
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intermediate  adjustment.  Consequently  the  rigidity  of  the  panel  would 
be  greatly  exaggerated.  The  decrease  in  the  effective  shear  modulus  of 
the  buckled  sheet  (reference  5)  because  of  the  larger  angle  subtended 
would  not  provide  a  sufficient  reduction  in  the  shear  rigidity.  The 
sheet  panel  was  therefore  taken  to  be  smaller,  the  natural  stopping  point 
being  the  position  of  the  stringer  next  to  the  most  highly  compressed 
stringer  in  the  actual  cylinder.  Thus  only  the  sheet  which  provides 
additional  stiffness  to  the  most  highly  compressed  stringer  is  considered. 
This  appears  consistent  with  the  assiimption  that  all  other  stringers  may 
be  neglected.  A  point  with  independent  freedom  of  motion  was  therefore 
considered  in  each  ring  at  the  intersection  with  the  free  edge  of  the 
sheet.  It  may  be  remarked  that,  if  the  rings  were  to  be  terminated 
there,  the  consequent  reduction  in  the  influence  coefficients  would  in 
general  be  negligibly  small. 

(it)  The  length  of  the  cylinder  was  considered  constant  and  equal  to 
six  times  the  ring  spacing.  For  the  three  cylinders  investigated,  PIBAL 
cylinder  10  and  GALCIT  cylinders  2$  and  65  (fig.  1),  this  corresponds  to 
1.5>  and  1.2  times  the  respective  diameters.  For  the  fuselage  of  a 

large  modern  transport  this  length  would  be  approximately  equal  to  the 
diameter.  The  following  table  may  be  set  up  on  the  basis  of  ejqperimental 
results  presented  in  the  references  given: 


Loading 

Limiting 

Increase  in  buckling  load 

value  of 
L'/B 

at  lower  value  of  L*/D 
(percent) 

Reference 

Compression 

1.5 

6  at  L'/D  =  1.0 

6 

Pure  bending 

2.0 

12  at  L'/D  =  1.2 

7 

The  buckling  load  is  practically  independent  of  the  length  if  the  length- 
to-diameter  ratio  L'/D  is  equal  to  or  larger  than  the  limiting  value 
given.  Examples  of  the  increase  in  the  buckling  loads  for  cylinders 
shorter  than  the  limiting  length  are  given  in  the  third  coliimn  of  the 
above  table.  The  length  assumed  in  the  calculation  will  be  in  general 
somewhat  shorter  than  the  limiting  length;  the  error  caused  by  this  may 
be  estimated  with  the  aid  of  the  above  table  to  be  at  most  10  or  l5  per¬ 
cent  of  the  buckling  load  of  a  cylinder  longer  than  the  limiting  length. 
The  effect  of  the  length  was  investigated  in  some  detail  with  test 
cylinder  25  of  the  GALCIT  series  (reference  8).  The  buckling  load  for 
this  cylinder  was  calculated  considering  different  numbers  of  bays  and 
the  results  are  shown  in  figure  i;.  It  may  be  seen  that  the  calculated 
buckling  loads  approach  some  constant  value  in  what  appears  to  be  an 
asymptotic  variation  and  that  the  difference  in  the  buckling  loads 
obtained  considering  six  or  eight  bays  is  small.  It  was  concluded  that 
the  small  improvement  in  accuracy  given  by  a  longer  structure  did  not 
warrant  the  increased  amount  of  work  required  to  obtain  it. 
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(5)  The  buckling  load  was  calculated  with  the  aid  of  the  simplified 
structure  for  test  cylinders  25  and  65  of  the  GALCIT  series  (reference  8) 
and  for  cylinder  10  of  the  PIBAL  series  (reference  9).  Those  specimens 
were  chosen  because  of  their  widely  different  characteristics  (see  fig.  1 
and  table  l).  Comparisons  of  the  results  of  the  present  analysis  with 
those  of  experiment  are  presented  in  table  IV  and  in  figures  k,  5,  and  6. 
The  calculated  buckling  loads  may  be  seen  to  be  consistently  higher  than 
the  corresponding  experimental  values.  The  percentage  errors  obtained 
are  not  considered  excessive,  however,  upon  comparison  with  the  results 
obtained  earlier  at  PIBAL  by  means  of  strain-energy  methods.  One  of  those 
solutions  (reference  lO)  gave  better  results  than  the  present  investi¬ 
gation,  but  required  a  prohibitive  amount  of  work. 

(6)  Approximate  deflected  shapes  .at  buckling  obtained  with  the  aid 

of  the  simplified  structure  of  figure  2  are  given  in  table  V  for  the  three 
cylinders  investigated.  The  same  table  also  gives  results  of  measurements 
made  on  some  actual  test  specimens  after  buckling  (reference  6).  It  may 
be  noticed  that  fair  agreement  has  been  obtained  between  measured  and 
calculated  values,  so  that  an  additional  indirect  experimental  check  has 
been  provided  on  the  reasonableness  of  the  simplified  structure  chosen. 

It  should  be  kept  in  mind  that  the  measurements  were  taken  after  the 
cylinders  had  buckled,  and  therefore  may  differ  from  the  actual  displace¬ 
ments  at  the  instant  of  buckling. 

(7)  Table  V  also  shows  that  the  radial  displacements  r  of  the  most 
highly  compressed  stringer  at  buckling  are  closely  represented  by 

r  =  sin^(nx/6L)  (It) 

where  n  =  ii,  5,  or  6.  The  rotations  may  be  closely  approximated  by 

m-5  =  (nTt/6L)  sin^"^(nx/6L)  cos  (nx/6L)  (5) 

The  values  obtained  with  n  =  5  represent  a  reasonable  average  of  all 
experimental  and  calculated  deflections,  and  therefore  this  value  of  n 
was  chosen  for  the  solution  with  assumed  displacements  which  was  described 
previously.  As  a  check,  the  buckling  load  of  GALCIT  cylinder  65  was  cal¬ 
culated  by  that  method  and  was  found  to  be  173Q  pounds,  the  buckling  load 
calculated  from  the  simplified  structure  without  the  assumption  of  dis¬ 
placements  was  1670  pounds  (fig.  5),  so  that  the  error  introduced  by  this 
assumption  is  only  3.6  percent  of  the  latter  value.  Cylinder  65  was 
chosen  since  it  is  the  specimen  for  which  the  agreement  between  assumed 
and  actual  displacements  is  the  poorest  (table  V).  It  should  be  remembered 
in  this  connection  that  it  was  shown  elsewhere  (reference  2)  that  the 
methods  of  calculation  used  in  this  report  are  not  too  sensitive  to  errors 
in  the  assumed  deflected  shape. 
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and  in  j  loSort”  " 

ivide  the  tenth  and  eleventh  rows  and  columns  by  L  and^  Process, 
rows  and  columns  (10)  through  (li,)  by  the  ^antlly  \ 

colu^  terMs  appearing  ip  the  lo»er  right-hand  corner'^.ay  Ln  he  »ritte'n 

some  function  of  kL.  The  ring  influence  coef¬ 
ficients  rr,  rm,  rt,  and  so  forth  are  equal  to  (El)  /d^  multinlied 

table  „in  he  a  /Jiji 

.  J>  ttl)atr„  1 


Gefftd^ 

(El)jX 


where  r  is  the  cylinder  radius  CftI  +hQ  k  j- 

in  its  0™  plane,  d  the  ciroun^erLl^  “ri^er  epaSnf^G  \h“" 

effective  shear  modulus  for  the  sheet  t  ^  ^ 

other  symbols  have  been  previously  dehned  The^eff  thickness,  and  the 
e^ensional  deformations  of  the  rings  J^necnLi  shearing  and 

the  two  additional  parameters-  ^  '  respectively,  are  represented  by 


I  =  A//Ar  I 

.  =  Ard7i^  J 

where  Ip  is  the  moment  of  inert  +Vi^ 

the  area  of  the  ring  and  A  *  •  tv,  cross  section,  Ap  is 

cross  section.  Reference  ll^'sho  ^  effective  shear  area  of  the  ring 
parameters  is  in  generS  nigUgSle.  ^^ese  two 

It  has  been  shown  in  reference  1?  that  th=  v  in- 
coque  cylinder  depends  on  the  para™,ter  A.  I„„  aSdiSnarpa^LetSr 
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were  established  in  that  the^Ss^highly  compressed  stringer 

is  the  buckling  strain  of  a  sheet  panel.  An  experi- 
at  failure  and  cj,  parameters  approximately 

mental  verification  of  the  f  ^  mnnncoaue  cylinders  is  given  in  refer- 
control  the  buckling  of  the  paraiiters  found  in  the  present 

ence  13-  It  may  be  seen  that  two  ^^ence  12,  while  the  param- 

development  are  the  same  since  the  shearing  rigidity  Ggff 

eter  T  includes  the  quantity  e/^cr’ 

was  found  in  reference  5  to  be  closely  approximated  by 


=  N  +  (1  -  N)e 
Gn 


-N  e/< 


where 


and  G, 


N  =  0.0275  [(2nr/d)  +  l] 
is  the  shear  modulus  of  the  sheet  material. 


curves  of  M.,  »hioh  cJ^lXrot 

the  parameters  A  insufficient  number  of  cylinders  is  _ 

reference  8  with  r/d  -  6.32.  curves  is  not  definitely  determned. 

available  so  that  the  posi  results  presented  do  not  contradict 

T+  mav  be  stated,  however,  that  tne  , 

lie  validity  of  the  four  parameters  established. 

CALCULATION  OF  BUCKLING  LOAD  OF  A  CYLINDER  WITH  A  CUTOUT 

Experimental  Investigation 

The  methods  developed  (piB^®cylindL^82^^fig.^8),  which 

consisted  of  a  thin  gg^gS^SLd  and  tested.  The  cutout 

evenly  spaced  rings  ^the  compression  side  of  the  cylinder, 

extended  circumferentially  for  ^  cylinder  through 

Pure  bending  moments  were  trans  .  test  rig  was  the  same  as 

heavy  rings  »l>lch  could  b®  Se^ncfu  This  cylinder  buckled 

that  used  in  the  cylinder  ^  ^®^gsged  stringer  (stringer  2  in  fig.  8) 

when  the  load  in  ‘‘"J®  J  ^gj.respoSded  to  a  total  applied  moment  of  ^ 

r58,S°iS:rJUs  pXUs  of  the  buckling  cylinder  are  shown  an 

figures  9  and  10. 
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Method  of  Calculation 

The  results  of  the  theoretical  investigation  indicated  that  the 
simplified  methods  evolved  for  complete  cylinders  are  not  satisfactory 
for  cylinders  with  cutouts.  For  such  cylinders,  the  distortion  under 
load  depends  primarily  on  the  geometry  of  the  cutout,  so  that  only  dis¬ 
placements  in  the  vicinity  of  the  cutout  require  consideration  in  the 
operations  table. 

With  consideration  of  the  symmetry  of  the  cylinder,  an  operations 
table  including  only  the  displacements  rgj_  and  r32  and  the  rotation 

mj.B2  (see  fig.  8)  was  set  up.  The  resulting  buckling  load  was,  for  all 

practical  purposes,  the  same  as  that  obtained  through  the  use  of  an 
operations  table  which  permitted  all  possible  generalized  displacements 
at  all  the  joints  of  the  cylinder.  Hence  the  assumption  that  only  the 
joints  in  the  vicinity  of  the  cutout  need  be  considered  in  the  operations 
table  was  justified  for  the  cylinder  with  a  cutout.  The  actual  experi¬ 
mental  deflected  shape  of  PIBAL  cylinder  82  (figs.  9  and  10)  shows  that 
the  major  distortions  took  place  in  the  vicinity  of  the  cutout,  and  that 
all  other  joints  may  be  assumed  to  have  had  zero  displacements. 

Cylinders  encountered  in  practice,  however,  will  be  of  a  more  compli. 
cated  construction  than  PIBAL  cylinder  82,  and  hence  the  simplified 
operations  table  described  above  may  not  be  sufficiently  complete. 
Depending  on  the  size  of  the  cutout,  it  is  suggested  that  the  operations 
table  be  expanded  so  as  to  include  all  the  joints  surrounding  the  par¬ 
ticular  cutout. 

The  buckling  load  obtained  for  PIBAL  cylinder  82  considering  only 
three  generalized  displacements  was  8I1.OO  pounds.  The  discrepancy  between 
the  theoretical  and  experimental  buckling  loads  was  attributed,  mainly, 
to  the  inaccuracy  of  the  value  of  the  effective  shear  modulus  Gg££  used 

in  the  calculations.  This  value  was  0.71Gq  and  was  taken  from  equa¬ 
tion  (8).  This  equation  is  based  on  bests  on  panels  buckled  because  of 
compression.  The  sheet  panels  in  the  present  cylinder,  however,  are 
imder  the  action  of  combined  compression  and  shear.  No  values  for  the 
effective  shear  modulus  of  curved  panels  under  such  a  loading  could  be 
found  in  the  literature,  but,  according  to  data  obtained  from  flat  panels 
(reference  lU),  it  appears  that  the  correct  value  of  Gg^.^.  should  be 

considerably  lower.  Furthermore,  as  is  shown  in  the  next  section,  there 
is  reason  to  reduce  the  shear  modulus  even  fiirther. 

The  calculations  for  PIBAL  cylinder  82  were  therefore  repeated  with 
an  assumed  value  of  Gg^^  =  O.IGq.  The  resulting  buckling  load  was 

5900  pounds,  which  may  be  seen  to  be  in  good  agreement  with  experiment. 
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Reduced  Effective  Shear  Modulus 

If  a  panel  of  sheet  is  not  in  a  buckled  state,  the  relation  between 
shear  stress  t  and  shear  strain  y  is  simply  Hooke's  law; 

x  =  Gq^  (9) 

If  the  panel  is  in  a  buckled  state,  a  relation  analogous  to  equation  (9) 
will  still  hold  between  the  average  shear  stress  and  the  average 

shear  strain  y^^,  provided  that  an  effective-  shear  modulus  Grgff  is 

used  in  place  of  Gq.  In  other  words. 


^eff^av 


and  Gg|>£  =  Gq  if  the  panel  is  not  buckled.  The  value  of  Gg|.^  wi 

represent  the  complex  state  of  stress  of  the  buckled  panel,  and  will 
presumably  vary  with  panel  dimensions  and  type  of  loading. 


will 


The  value  of  Gg^^  is  the  proportionality  factor  between  the  average 

shear  stress  and  the  average  shear  strain.  In  problems  of  instability, 
however,  it  is  desired  to  know  the  relation  between  a  small  increase  in 
stress  and  a  small  increase  in  strain  d.(7ay)-  This  relation 

will  again  have  the  same  form  as  equation  (9),  if  only  a  reduced  effective 
shear  modulus  *Ieffj.ed  used  in  place  of  Gq.  In  other  words, 

‘^('^av)  ~  ^effped  ^("^a-v) 

Thus  this  new  modulus  represents  the  resistance  the  panel  will  offer 
against  distortions  additional  to  those  represented  by  y^^.  According 

to  the  previous  discussion,  this  new  modulus  will  also  depend  upon  the 
dimensions  of  the  panel  and  upon  the  amount  of  shearing  and  compressive 
loads  present. 

If  equation  (10)  is  written  in  differential  form  as 


d(Gpf.f) 

('^av)  =  ^ 

V  av; 


comparison  with  equation  (11 )  indicates  that 
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The  following  remarks  may  be  made  about  the  reduced  effective  shear 
modulus 

(1)  Geff^g^  =  Gq  when  the  panel  is  not  in  a  buckled  state 

(2)  ^effpg(^  ^  ^eff  when  the  average  shearing  strain  in  the  panel 
is  zero  immediately  before  buckling 

(3)  ‘^effj;.ed  “  *^eff  when  there  is  no  change  of  shearing  strain 
during  buckling  of  the  structure  \inder  consideration 

(li)  Geffj^gi^  <  Ggff  in  all  other  cases,  since  in  general  the  modulus 
Geff  decreases  with  increasing  shear  strain  so  that  the  second 

term  in  the  right-hand  side  of  equation  (13)  is  negative 

The  latter  case  applies  to  the  cylinder  with  a  cutout.  The  low 
value  assumed  for  the  shear  modulus  in  the  calculations  is  therefore 
plausible . 


Polytechnic  Institute  of  Brooklyn 

Brooklyn,  N.  Y.,  August  31,  19i|.8 
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APPENDIX  A 

BASIC  THEORI 


The  procedures  developed  in  this  report  for  the  calculation  of  the 
buckling  load  of  a  monocoque  cylinder  are  based  on  methods  developed  in 
reference  2  which  make  use  of  an  operations  table  similar  to  that  used 
in  Southwell's  method  of  systematic  relaxation.  These  methods  are  out¬ 
lined  here,  more  rigorous  proofs  being  given  in  reference  2.  Rigorous 
proofs  are  only  given  here  for  some  modifications  of  these  methods  which 
were  not  discussed  in  that  reference. 

Consider  several  points  in  the  structure  in  question  distributed  so 
as  to  cover  the  entire  structure.  Let  these  points  be  numbered  consecu¬ 
tively  from  1  to  n.  The  generalized  force  exerted  on  joint  i  by  a  gener¬ 
alized  displacement  xj  at  joint  j  (all  joints  but  j  being  considered 
temporarily  rigidly  fixed)  may  be  denoted  by  apjXj.  The  quantity  aj^j 
is  called  an  influence  coefficient.  If  Fj_  is  the  generalized  external 
force  acting  at  joint  i,  the  equilibrium  condition  for  the  ith  joint  is 

n 

(Al) 

0-1 

provided  that  the  principle  of  superposition  is  valid.  If  equation  (Al) 
is  written  for  every  joint  in  the  structure,  a  set  of  linear  simultaneous 
equations  will  result  with  generalized  displacements  as  unknowns.  The 
array,  or  matrix,  of  the  coefficients  of  this  set  of  equations  is  called 


the  operations  table  and  may  be 

written 

as 

^11  ^12  •  • 

.  a]Li  • 

aij  .  . 

•  ®-ln 

o 

a2i  a.22  •  • 

•  ^21  ■ 

•  • 

«  # 

a2j  .  . 

•  a2n 

•  •  • 

= 

^il  ^i2  •  • 

.  aii  . 

•  • 

•  • 

aij  .  . 

•  •  • 

•  ^In 

ajp  aj2  •  • 

.  a^i  . 

•  • 

•  • 

•  ^jn 

^n2  •  • 

•  ^ni  • 

•  • 

a^j  .  . 

•  ^nn 

As  a  consequence  of  Maxwell's  reciprocal  theorem  ^ij  =aji- 
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Equation  (Al)  represents  the  equilibrium  conditions  for  the  given 
structure  in  terms  of  displacements.  In  general^  the  determinant  Aj^ 

will  not  be  equal  to  zeroj  then^  only  one  set  of  displacements  may  be 
found  which  will  satisfy  the  equilibrium  conditions.  If  the  determinant 
Aj^^  vanishes,  however,  more  than  one  such  set  of  displacements  will  exist. 
This  is  physically  possible  only  at  neutral  equilibrium,  or,  which  is 
the  same,  at  a  buckling  load.  This  leads  to  what  was  called  in  refer¬ 
ence  2  the  Determinant  Method,  the  basis  of  which  is  the  fact  that  the 
lowest  load  at  which  the  determinant  A„  vanishes  is  the  lowest  buckline 
load.  ^ 


A  proof  will  now  be  given  of  the  fact  that  in  general  at  the  lowest 
buckling  load  the  last  main— diagonal  element  of  the  auxiliary  matrix  of 
the  method  for  evaluating  determinants  given  in  reference  3  is  equal  to 
zero.  Let  the  symbol  Aj_  stand  for  the  determinant 


ail  ai2  •  .  . 

a2i  a22  .  .  .  a2i 


an  ai2  •  •  • 


(A3) 


The  value  of  this  determinant  is  equal  to  the  product  of  the  first  i  main 
diagonal  terms  of  the  auxiliary  matrix.  If  auxiliary— matrix  elements  are 
denoted  by  the  symbol  aqj,  then 


Ai  -  ajj  -  “11  “22 

J=1 


“kk  •  *  •  “ii 


(Aii) 


Let  Aj^  be  the  first  of  these  determinants  to  vanish]  then  theorem  2 
of  reference  2  gives 


Ak  =  Ak+i  =  .  .  .  =  An  =  0  (A3) 

where  An  is  the  determinant  given  in  equation  (A2).  Two  cases  may 
then  be  considered; 

Case  1;  k  =  n.-  In  the  case  where  k  =  n,  Aj^  is  the  only  one  of 
these  determinants  which  vanishes.  By  equation  (aU)  the  only  factor 
which  is  contained  in  A^  and  in  no  other  Aj_  determinant  is  Onn? 
which,  as  was  to  be  proved,  must  therefore  vanish. 

Case  2;  k  <  n.-  In  the  case  where  k  <  n,  the  method  of  reference  3 
fails  to  give  any  terms  beyond  which  of  course  is  zero.  Here  a 
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is  obviously  not  the  first  term  to  vanish,  but  in  this  case  several  of 
the  higher  buckling  loads  are  identical  with  the  first.  This  case  is 
expected  to  occur  rather  infrecjuently. 

The  Energy  Method' of  reference  2  is  based  on  the  condition  that  the 
second  variation  of  the  total  potential  energy  must  vanish  at  buckling. 
This  condition  may  be  written  as 


Q 


1=1  j=l 


aijXiXj  =  0 


(A6) 


This  equation  is  satisfied  at  buckling  by  the  buckling  displacements. 
In  the  Energy  Method  some  of  these  displacements,  say  Xp,  Xp+]_ 


"nj 


are  guessed;  then  the  others  are  obtained  from  the  conditions 


=0  k  =  1,  2,  .  .  .,  p  -  1  (A7) 

dXk 


The  matrix  of  the  coefficients  of  these  simultaneous  equations,  including 
the  constant  terms,  is  the  reduced  matrix  A',  where: 


A' 


ail 

• 

C\J 

1 — 1 
cti 

^21 

^22 

•  •  • 

ap-1,1 

ap-1,2  • 

a'p,l 

^'p,2.  • 

ai,p-l  ^'l,p 

^2,p-l  ^'2,p 


a-p-l,p-l  a'p-i,p 
a'p,p-l  a'p^p 


(a8) 


in  which 


iP 


n 


J=P 


and 


^'p,p 


n  n 

=  ZT  ZI 

i=P  D=P 


It  will  now  be  proved  that  the  Determinant  Method  may  be  applied  to  the 
operations  table  of  equation  (A8 ) . 
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It  is  well-known  (see  reference  l5)  that  any  quadratic  form  Q  may  be 
put  in  the  form 


~  zZ  (A9) 

i=l 

where  the  bj_  quantities  are  constants,  and 


n 


where  the  cj  quantities  are  constants.  It  is  assumed  that  there  are 

n  linearly  independent  quantities  L^.  This  assumption  entails  no  loss 

of  generality  since  in  the  case  in  which  it  is  not  true  some  of  the 
constants  b^  will  be  zero. 


^y  means  of  equation  (A9)  and  table  2  of  reference  2  the  followine 
table  may  be  set  up: 


Sign  of  bq 

Sign  of  Q 

Classification  of  Q 

Type  of 
equilibrium 

All^  less 
than  zero 

Q  <  0  always 

(Q  -  0  if  xj  *  0) 

Negative  definite 
nonsingular 

Stable 

Varying;  none, 
zero 

Q  may  be  positive, 
negative,  or  zero 

Indefinite 

nonsingular 

Unstable 

All,  greater 
than  zero 

Q  >  0  always 

(Q  .  0  If  xj  «  0) 

Positive  definite 
nonsingular 

Unstable 

Some,  zero;  all 
others,  negative 

Q  =  0  or 

Q  <  0 

Negative  definite 
singular 

Neutral 

Some,  zero; 
others,  vaiying 

Q  may  be  positive, 
negative,  or  zero 

Indefinite  singular 

Unstable 

Some,  zero;  all 
others,  positive 

Q  =  0  or 

Q  >  0 

Positive  definite 
singular 

Unstable 

_ i 

% 
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Two  conditions  for  neutral  equilibrium  have  been  thus  set  up: 

(1)  The  vanishing  of  the  determinant  Aj^  of  the  quadratic  form  Q 

(2)  The  fact  that  some  of  the  b^'s  equal  zero,  while  all  others 
are  negative 

As  both  conditions  are  necessary  and  sufficient,  they  are  equivalent  and 
may  be  used  interchangeably. 

When  some  of  the  displacements  are  guessed  as  previously  explained, 
the  quadratic  form  Q  becomes  the  reduced  quadratic  form  Q',  where 

Q' 


where  the  by'  quantities  are 


Li' 


constants  and 
P-1 

=  Cp'  +  ^  Cj'Xj 


(All) 


(A12) 


The  statements  of  the  above  table  may  now  be  applied  to  the^  quadratic 
form  since  is  the  value  the  quadratic  form  Q  will  take  on 

when  the  displacements  Xp,  .  .  .  x^  are  assumed.  Neutral  equilibrium 
will  then  exist  when  some  of  the  bj^  *  s  are  zero  and  all  others  are  neg¬ 
ative.  But  this  condition  is  equivalent  to  the  condition  that  the  determi¬ 
nant  corresponding  to  the  quadratic  form  Q*  must  vanish.  This  determinant 
is  obtained  by  multiplying  out  the  right-hand  side  of  equation  (All)  and 
expressing  the  x'esult  in  the  form 


p-1  p-1  p-1  P;:! 

Q'  =  aqjXiXj  a'ipXq  +  ^  a'pqxq  +  a'pp  (A13) 

i=l  j=l  i=l 

The  determinant  corresponding  to  Q'  will  then  be  seen  to  be  identical 
with  A'  of  equation  (A8). 

It  therefore  follows  that  the  vanishing  of  the  determinant  A'  of 
equation  (A8)  corresponds  to  neutral  equilibrium.  If  the  displacements 
X  ,  •  •  'j  were  not  chosen  exactly  equal  to  the  displacements  of 

the  structtire  at  buckling,  an  approximate  value  of  the  buckling  load  will 
be  obtained  rather  than  an  exact  one.  It  was  proved  in  reference  2  that 
this  approximate  load  will  be  higher  than  the  actual  one. 
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APPENDIX  B 

GENERAL  FORMULAS  AND  NUMERICAL  EXAMPLE 


In  this  appendix  is  presented  the  procedure  for  the  determination 
of  the  influence  coefficients  required  in  setting  up  the  operations 
tables  II  and  III.  Since  many  of  the  formulas  used  in  the  analysis  are 
scattered  throughout  the  literature,  some  of  these  are  given  here, 
together  with  appropriate  reference.  When  a  formula  is  not  listed,  refer¬ 
ence  to  its  source  is  given.  A  numerical  example  is  also  given  illus¬ 
trating  both  methods  suggested  for  the  calculation  of  the  buckling  loads 
of  cylinders  without  cutout. 


Influence  Coefficients 

The  operations  tables  (tables  II  and  III)  contain  three  types  of 
influence  coefficients, .which  represent  the  effects  of  the  rings,  the 
sheet  covering,  and  the  stringers.  The  methods  by  which  each  type  is 
determined  are  outlined  below. 

Ring  influence  coefficients.-  The  ring  influence  coefficients  are 
characterized  by  the  symbol  as,  for  example,  rrj/[,  rtj^,  or  tnp. 

These  coefficients  may  be  determined  from  reference  11,  in  which  they 
may  be  seen  to  depend  on  the  three  parameters  p,  y,  and  5,  where  p 
is  the  central  angle  of  the  ring  segment,  and  y  and  ^  are  defined  in 
equations  (7).  Any  one  of  the  three  following  ways  may  be  used  for  the 
calculation  of  these  coefficients: 


(l)  General  formulas  are  given  in  equations  (20),  (27),  (28),  and 
(29)  of  reference  11.  As  an  example,  the  formula  for  rtpj  is  repeated 
here  with  a  slight  change  in  notation: 


-rt 


d3 


W  (El): 


p  d3 

-‘“iiTirj 


_  _  n3  ^tr  r-^ 
r  ^  A  (El)r 


(Bl) 


-^tr  +  (l/r)  [(1/1)  -  f|ftr2 


where 
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and 

=  fAl  +  (1/r)  [fA2  +  (l/l)(fA3)]  ^  (1/t)  +  [(1/^)  -  ^^fAU 

in  which  the  quantities  f  are  functions  of  the  central  angle  p  and 
are  given  in  reference  11  by  equations  (2k)  and,  for  specific  values 
of  p,  by  table  II  and  figures  3  to  13- 

(2)  Influence  coefficients  obtained  from  the  formulas  mentioned 
under  item  (l)  are  tabulated  in  tables  III  and  IV  of  reference  11  for 
specific  values  of  the  parameters  p,  y,  and  I. 

(3)  From  the  tabulated  values  mentioned  in  item  (2)  above,  curves 
were  plotted  which  are  presented  in  figures  lli  to  8$  of  reference  11. 

It  should  be  noted  that  in  general  these  curves  are  accurate  only  for 
values  of  p  =  l5° • 

Sheet  influence  coefficients.-  The  shear  in  the  sheet  covering  is 
represented  by  the  influence  coefficients  containing  the  quantities  a 
and  Aj,  as,  for  example,  2aj.^Aj  or  2a^ajjAj  dj.  The  quantities  a 

may  be  determined  in  good  approximation  from  the  formulas 
Op  =  O.lp 

an,  =  -0.5(1  -  0.01666  •  •  •  p2) 

Un  =  -0.008333  ’  •  •  p(l  0.011i286p^) 


"bakGri  from  pcigo  27  of  reforGncG  11.  It  should  bG  notGd  "bhat  in  this 


rof  GPGncG 


and  Ur,  arG  dGnotGd  as 


nq/(  L^q ),  rGspGctivGly . 

ThG  quantity  Aj  is  givGn  by 


Aj  .  Oeff  ^  (B3) 

where  the  effective  shear  modulus  Ggff  is  given  by  equation  (8)  of  the 
prGSGnt  rGport. 
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The  buckling  strain  of  a  sheet  panel  appearing  in  this  equa¬ 

tion  was  obtained  by  means  of  Redshaw's  formula: 


(eflat/2^  +  \(€flat/2)^  + 


curved 


where 


«flat 


12  1  -  v 


as  shorn  in  reference  16,  in  which  k'  is  the  end-fixity  coefficient 
V  IS  Poisson's  ratio,  and  ^  xxxcxexio. 


Wved  ■  0-6(t/r)  (B6) 

1  +  0.00i;(E/Fcy)  ^ 

as  shown  in  reference  17,  in  which  is  the  yield  stress  of  the 

material. 

The  effect  of  normal  stresses  in  the  sheet  covering  is  taken  into 
account  by  an  effective  width  of  sheet  as  discussed  in  the  next  section. 

Stringer  beam-column  influence  coefficients.-  The  beam-column 
effectsin  the  stringers  are  represented  by  the  influence  coefficients 
containing  the  load  P,  as,  for  example,  (p/Di)(1  -  c)  or  Pks/Di 

f  of  beam-column  influence  coefficients  is  presented  in 

tables  VI,  VII,  and  VIII  in  which  the  sign  convention  as  well  as  defi- 
nitions  of  symbols  are  given.  All  these  coefficients  are  functions  of 
the  quantity ^  k  defined  in  equation  (2).  The  following  formula  which 
was  derived  in  reference  18  on  the  basis  of  work  contained  in  reference  IQ 
IS  suggested  for  the  calculation  of  the  effective  width  2w: 

2w  =  (lA)(d/r)  jo.lt  +  1.535  [(t/d)(r6  -  0.3t)rV2]^/^| 

If  the  load  P  causes  a  stringer  stress  which  is  higher  than  the 
proportional  limit  of  the  material,  the  modulus  of  elasticity  £  which 
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appears  in  equation  (2)  must  be  reduced  in  an  appropriate  manner 
this  investigation  Von  KannAn's  formula  was  used; 

EEt/(\lE  + 

where  is  the  tangent  modTilus. 


Numerical  Example 

As  an  example  of  the  application  of  the  methods  sugpsted  for  the 
calculation  of  the  buckling  load  of  cylinders  without  cutout  the  buckling 
load  is  determined  here  for  GALCIT  cylinder  65  of  reference  8.  The 
characteristics  of  this  cylinder  are  given  in  figure  1.  Table  I  is 
readily  set  up  with  the  aid  of  these  characteristics  and  the  equations 
SSed  in  the  previous  section.  The  influence  coefficients  required 
for  tables  II  or  III  can  then  be  calculated. 

Ring  influence  coefficients.-  The  values  of  the  ring  influence 
coefficients  corresponding  to  pj,  yj,  and  Ij  (see  table  l)  happen 
to  appear  la  tablee  III  of  reference  11.  Howe^r, 

do  not  contain  coefficients  corresponding  to  ijj,  it  is  necessary 

use  values  interpolated  from  the  appropriate  curves  of  reference  11. 

The  following  values  were  obtained  for  the  ring  influence  coefficients: 


Ring  segment  I 


6.393 


Ring  segment  II 


MiVEI 


=  -20.38 


=  -33.0 


%\EI 


%l\  El 


‘MiiVk 


II\  Elr 


=  170 
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The  beam-column  influence  coefficients  required  for  table  II  are 

then 


=  281.31 

1^1 


~  .-g^  =  1337.6 
°1 

No  beam-column  influence  coefficients  need  be  calculated  if  table  III 
and  figure  3  are  used. 

Calculation  of  buckling  load  by  table  II.—  The  state  of  stability 
of  the  structure  at  a  load  of  F^T^^o’pounds  may  now  be  investigated 
by  introducing  all  the  above  influence  coefficients  into  table  II  and 
evaluating  the  corresponding  determinant.  The  results  corresponding  to 
the  assumed  load  of  1^50  pounds,  as  well  as  those  for  the  loads  of  I650, 
1680,  and  1713  pounds,  are  presented  in  figure  3*  In  this  figure  are 
plotted  both  the  values  of  the  determinant  and  the  values  of  the  last 
main-diagonal  term  Unn  of  the  auxiliary  matrix  of  reference  3.  The 
intercept  of  the  curves  in  this  figure  may  be  read  off  and  corresponds 
to  the  buckling  load.  It  should  be  noted  that  only  the  stringer  beam- 
column  influence  coefficients  vary  as  the  assumed  load  is  changed. 

Calculation  of  buckling  load  by  table  III.—  The  value  of  the 
quantity  K  of  equation  (3)  was  found  by  the  method  of  reference  3  to 
be  -3.ii72.  From  figure  3,  this  value  may  be  seen  to  correspond  to 
kXj  =  from  which  =  1730  pounds.  Care  must  be  taken  that  a 

reduced  modulus  (equation  (B8))  be  used  if  the  stress  at  buckling  is 
above  the  proportional  limit. 


P(hL  -  .s)_  ^  223U.6 
kDi 

^  (s  -  kLc)  =  6191.7 
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APPENDIX  C 

SIMBOLS 

determinants  or  matrices 

effective  cross-sectional  area  of  stringer 
effective  cross-sectional  area  of  ring 
effective  shear  area  of  ring  cross  section 

reduced  matrix 
rings 

Young's  modulus 

bending  rigidity  of  a  ring  in  its  ovm  plane 

radial  bending  rigidity  of  stringer  and  its 
effective  width  of  sheet 

reduced  modulus 
tangent  modulus 
yield-point  stress 
generalized  external  force 
effective  shear  modulus 
reduced  effective  shear  modulus 

shear  modulus  of  sheet  material  under  no  com¬ 
pressive  load 

moment  of  Inertia  of  ring  cross  section  plus 
effective  width 


K 

L 

Li 

L'/D 

Mt 

N 


P 

Per 


Q 

Q' 

R 

T 


s' ip 
bi,  bi' 

d 

f 

f(kL) 
i>  J 


moment  of  inertia  of  stringer  plus  effective  width  of 
sheet  for  radial  bending 

constant 

ring  spacing 

linear  function 

ratio  of  total  cylinder  length  to  cylinder  diameter 

moment  causing  bending  of  stringer  (vector  pointing 
in  tangential  direction) 

moment  causing  bending  of  ring  in  its  plane;  also 

0.0275  [(2Tir/d)  +  ij 

axial  stringer  load 

load  in  the  most  highly  compressed  stringer  at  the  instant 
of  buckling 

quadratic  form 

reduced  quadratic  form 

radial  force 

tangential  force 

element  of  operations  table 

element  of  last  row  or  column  of  reduced  matrix 

constants 

constants 

circumferential  stringer  spacing 
function  of  p 
function  of  kL 
indices . 


k  =\|P/(El)strr 
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end-fixity  coefficient 

rotation  causing  bending  of  stringer  (vector  pointing 
in  tangential  direction) 

niimber  of  stringers;  number  of  generalized  displacements; 
rotation  causing  bending  of  ring  in  its  plane 

index 

radius;  radial  displacement 
ring  influence  coefficients 

sheet  thickness;  tangential  displacement 
effective  width  of  sheet 
longitudinal  axis  of  stringer 
generalized  displacement 
parameter 
parameter 

parameter 

auxiliary-matrix  element 

last  main-diagonal  element  of  auxiliary  matrix 

functions  of  p  required  for  sheet  influence 
coefficients 

central  angle  of  a  ring  segment  (d/r) 
parameter;  shear  strain 

strain  in  most  highly  compressed  stringer  at  failure 
buckling  strain  of  a  sheet  panel 

buckling  strain  of  nonreinforced  circular  cylinder 
under  uniform  axial  compression 


26 
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^flat 

V 

I 

a 

T 

Subscripts: 

A,  B,  C,  D 

F 

M 

av 

exp 

1,  2,  3 
I,  II 


buckling  strain  of  flat  panel  under  uniform  compression 

Poisson *s  ratio 
parameter 

compressive  stringer  stress 
shear  stress 

rings 

fixed 

movable 

average 

experiment 

stringers 

regions 
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All  cylinders  are  aluminum  (E  =  10. 5  x  10^  psi,  =  3.9  x  lo^  psi 
V  =  0.3,  E-t  from  reference  20).  ^ 

2 

In  general  €  must  be  guessed.  In  the  present  investigation  ^  Wi 
available  from  experimental  results  (references  8  and  9). 

not  required  since  ^fiat  ^curved  given  in  reference 


TABLE  n 

OPERATIONS  TABLE  FOR  THE  SIMPLIFIED  STRUCTURE 
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NACA 


TABLE  HI 

OPERATIONS  TABLE  FOR  SOLUTION  WITH  ASSUMED  DISPLACEMENTS 
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TABLE  IV 


COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL  BUCKLING  LOADS 


Cylinder 

Method  of 
calculation 

Calculated 

^cr 

(lb) 

Experimental 

^cr 

(lb) 

Difference 

(percent) 

GALCIT  25 

Table  II 

955 

769 

2li.2 

GAICIT  65 

Table  II 

1670 

1371 

21.8 

GALCIT  65 

Table  III 

1730 

1371 

26.2 

PIBAl  10 

Table  II 

1*850 

3751* 

29.2 

TABLE  V 

ASSUMED,  CALCULATED,  AND  EXPERIMENTAL  DEFLECTED  SHAPES 
[For  sign  convention  and  nomenclature  see  fig.  2] 


Shape 

^B 

^D 

mg 

JHq 

0.06250 

0.56250 

1 

0.22672(1/L) 

0.680l6(lA) 

0.031250 

0.1»8713 

1 

0. 11*170  (lA) 

0. 73636  (lA) 

0.015625 

0.1*2188 

1 

0. 085020  (lA) 

0.76520(1  A) 

PIBAL  oyl.  10 
(Calculated) 

-0.01691*9 

0.50229 

1 

0. 15086  (lA) 

0. 80011  (lA) 

GALCIT  cyl.  25 
(Calculated) 

0.066323 

0.1*8355 

1 

0.093608(1  A) 

0.82808  (lA) 

GALCIT  cyl.  65 

(Calculated) 

-0,16089 

0.3061*0 

1 

-0,21075(1A) 

1.1631(1A) 

GALCIT  cyl.  25 
(Experimental) 

0.222 

0.667 

1 

- ' — 

GALCIT  cyl.  27 
(Experimental ) 

0.0370 

■ 

0.201* 

1 

- - 

— 

GALCIT  cyl.  30 
(Experimental) 

0.0588 

■ 

o.ia2 

• 

1 

— 

- - - 

GALCIT  cyl.  35 
(Experimental) 

0.0893 

0. 1*1*6 

1 
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2(1  -  c )  -  kLs 

s 

=  sin  kL 

1  -  (kL/t) 

c 

=  cos  kL 

p/ei 

t 

=  tan  kL 

Sign  convention 

. 

Forces  on  constraints 

positive  as  shown 


3h 
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TABLE  VII 

BEAM-COLUMN  INFLUENCE  COEFFICIENTS 
[Axial  end  load  P  equal  to  zero”] 


Forces 

Displacement 

Unit 

Fa 

^A 

Fb 

Mb 

6a 

Both  ends 
fixed 

-12^ 

l3 

A  El 

12  FI 
l3 

A  El 

■"15 

“■A 

-6^ 

l2 

-1,11 

L 

-2  ^ 

L 

12 

L^ 

A  El 

A  El 

"? 

ms 

A  El 

L 

A  El 

-^T 

6a  . 

End  A 
fixed 

-3^ 

L^ 

-3^ 

l2 

3FI 

l3 

0 

“A 

End  B 
pin- 
jointed 

-3  ^ 

^  2 

0 

6b 

-3^ 

l3 

0 

6a 

Both 

ends 

pin- 

jointed 

0 

0 

0 

0 

6b 

0 

0 

0 

0 

Sign  convention 
Forces  on  constraints 


X 


/ 
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TABLE  VIII 

BEAM-COLUMN  INFLUENCE  COEFFICIENTS 
[Tensior^ 


Displacement 

Unit 

Forces 

(1) 

Fa 

”a 

Fb 

^B 

6a 

Both  ends 
fixed 

p  ks' 

Di' 

^l' 

P  ks' 

Dl' 

(C  -  Ij 

Di' 

®A 

P  (C  -  1) 
Di' 

mssi 

p  (c'  -  1) 
Di' 

MB 

pkil 

p  (C  -  1) 

Dl' 

p  ks' 

Dl’ 

p  (C  -  1) 

Dl* 

®B 

P  (C  -  1) 
Di' 

P  (c  -  1) 

6a 

End  A 

P  ^ 

-p  1 

P  . k 

0 

fixed 

^  t'D2’ 

D2‘ 

"  t’D2’ 

% 

End  B 

-p  — 

D2' 

-p-i. 

D2' 

0 

pin- 

jointed 

P  ^ 

p 

P  k  _ 

0 

^  t’D2’ 

^  Dg' 

‘  t'Dp' 

6a 

Both 

ends 

__  P 

"l 

0 

F 

L 

0 

6b 

pin- 

jointed 

P 

L 

0 

P 

^  L 

0 

1di '  =  2(1  -  c ' )  +  kLs ' 
Dg'  =  (kL/t)  -  1  ‘ 
k2  =  P/EI 


s’  =  sinh  kL 
c ’  =  cosh  kL 
t’  =  tanh  kL 


Sign  convention 
Forces  on  constraints 


P,  positive  as  shovn 


PIBAL  cylinder  10 


t  =  0.010' 
^  17S-T 


In 


Figure  2.-  Simplified  structure  for  setup  of  operations  tables. 
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Figure  4.-  Determination  of  buckling  load  of  GALCIT  cylinder  25  for  two,  four 

six,  and  eight  bays. 


10  X  loS 


P,  lb 


Figure  5.-  Determination  of  buckling  load  of  GALCIT  cylinder  65  by  table  II, 


P,  lb 


Figure  6.-  Determination  of  buckling  load  of  PIBAL  cylinder  10  by  table  II, 
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Figure  9.- 


Side  view  of  PIBAL  cylinder  82  after  buckling. 


